CYCLIC CLASSES AND ATTRACTION CONES IN MAX ALGEBRA 



SERGEI SERGEEV 

Abstract. In max algebra it is well-known that the sequence A'^, with A an irreducible 
square matrix, becomes periodic at sufficiently large k. This raises a number of questions 
on the periodic regime of A'' and A'' (g) x, for a given vector x. Also, this leads to the 
concept of attraction cones in max algebra, by which we mean sets of vectors with 
ultimate orbit period not exceeding a given number. 

This paper shows that some of these questions can be solved by matrix squaring 
{AjA^jA^, ...), analogously to recent findings of Semancfkova [39l [40] concerning the 
orbit period in max-min algebra. Hence the computational complexity of such problems 
is O(n^logn). The main idea is to apply an appropriate diagonal similarity scaling 
A ^ X^^AX, called visualization scaling, and to study the role of cyclic classes of the 
critical graph. 

For powers of a visualized matrix in the periodic regime, we observe remarkable sym- 
metry described by circulants and their rectangular generalizations. We exploit this 
symmetry to derive a system of equations for attraction cone, and present an algorithm 
which computes the coefficients of the system. 



1. Introduction 

By max algebra we understand the analogue of linear algebra developed over the max- 
times semiring Mmax,x which is the set of nonnegative numbers M+ equipped with the 
operations of "addition" a©6 := max(a, b) and the ordinary multiplication a®b := axb. 
Zero and unity of this semiring coincide with the usual and 1. The operations of the 
semiring are extended to the nonnegative matrices and vectors in the same way as in 
conventional linear algebra. That is if A = (ajj), B = (bij) and C = ( matrices of 

compatible sizes with entries from M+, we write C = A (B B ii Cij = aij © bij for all i,j 
and C = A0 B if Cij = ^® aikbkj = maxkiO'ikbkj) for all If A is a square matrix over 
M+ then the iterated product A^ A^ ...^ A in which the symbol A appears k times will 
be denoted by A'^. 
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The max-plus semiring Minax,+ = (K U {— oo},© = max,® = +), developed over the 
set of real numbers M with adjoined element —oo and the ordinary addition playing the 
role of multiplication, is another isomorphic "realization" of max algebra. In particular, 
X I— > exp(a;) yields an isomorphism between ]Rmax,+ and Mmax.x- In the max-plus setting, 
the zero element is — oo and the unity is 0. 

The min-plus semiring Mmin,+ = (IR U {+00},© = min, © = +) is also isomorphic to 
Rmax,+ and Mmax.x- Another well-known semiring is the max-min semiring Mmax,min = 
(M U {—00} U {00}, © = max, © = min), see [25l [39l HQ], but it is not isomorphic to any 
of the semirings above. 

Max algebraic column spans of nonnegative matrices A G M"^" are sets of max linear 
combinations of columns ©"^^ otiA.i with nonnegative coefficients Oj. Such column spans 
are max cones, meaning that they are closed under componentwise maximum © and 
multiplication by nonnegative scalars. There are important analogies and links between 
max cones and convex cones p3l [T7| W2[ W\\ . 

The maximum cycle geometric mean X{A), see below for exact definition, is one of the 
most important charasteristics of a matrix A G M"^" in max algebra. In particular, it 
is the largest eigenvalue of the spectral problem A © x = Ax. The cycles at which this 
maximum geometric mean is attained, are called critical. Further, one consideres the 
critical graph C{A) which consists of all nodes and edges that belong to the critical cycles. 
This graph is crucial for the description of eigenvectors [HI [151 EE] . 

The well-known cyclicity theorem states that if A is irreducible, then the sequence A^ 
becomes periodic after some finite transient time, and that the ultimate period of A^ is 
equal to the cyclicity of the critical graph [3l |15l [26]. Generahzations to reducible case, 
computational complexity issues and important special cases of this result have been 
extensively studied in [ISl 123 [251 131 [S3] • 

In this paper we study the behaviour of matrix powers and orbits A'^ © x in the ir- 
reducible case in the periodic regime, i.e., after the periodicity is reached. One of the 
main ideas is to study the periodicity of visualized matrices, meaning matrices with all 
entries less than or equal to the maximum cycle geometric mean. This study provides a 
connection to the theory of Boolean matrices [6|, [29] . 

In Boolean matrix algebra, one considers components of imprimitivity of a matrix 
[ni[2n], or equivalently, cyclic classes of the associated digraph [1]. In max algebra, cyclic 
classes of the critical graph have been considered as an important tool in the proof of the 
cyclicity theorem mentioned above, see [26] Sect. 3.1. Recently, the cyclic classes appeared 
in max-min algebra [39l HQ], where they were used to study the ultimate periods of orbits 
and other periodicity problems. It was shown that such questions can be solved by matrix 
squaring {A, A'^, A"^, A^, ...), which yields computational complexity 0(n'^ log n). 
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We show that the problems of computing ultimate period and matrix powers in the 
periodic regime can be solved by matrix squaring in max algebra, which yields the same 
complexity bound O(n^logn). This is achieved by exploiting visualization, and cyclic 
classes of the critical graph. Further it turns out that the periodic powers of visualized 
matrices have remarkable symmetry described by circulant matrices and their rectangular 
generalizations. We use this symmetry to derive a system of equations for attraction cone, 
meaning the max cone which consists of all vectors x whose ultimate period of A^^x does 
not exceed a given number. We also describe extremals of attraction cones and present an 
algorithm for computing the coefficients of this system in the case when C{A) is strongly 
connected. 

The contents of the paper are as follows. In Section 2 we revise two important topics 
in max algebra, namely the spectral problem and Kleene stars. In Section 3, we speak of 
the visualization and the connection to the theory of Boolean matrices which it provides, 
see Propositions 13.11 and 13. 3[ In Section 4, we study basic properties of matrix powers 
in the periodic regime, see Propositions 14.51 - 14. 7[ The problems which can be solved by 
matrix squaring are described in Theorem l4.11[ In Section 5 we introduce some useful con- 
structions associated with irreducible visualized matrices and their powers, namely, core 
matrix, CS'-R-representation, and describe their circulant simmetries in Proposition 15.41 
In Section 6 we derive a concise system for attraction cone, see Theorem 16.21 We also 
describe extremals and present an algorithm for computing the coefficients of this system 
in the case when C{A) is strongly connected. We conclude with Section 7 which is devoted 
to numerical examples. 

As Minax,+ and Mmax,x are isomorphic, we use the possibility to switch between them, 
but only when it is really convenient. Thus, while the theoretical results are obtained over 
max-times semiring, which looks more natural in connection with diagonal matrix scaling 
and Boolean matrices, the examples in Section 7 are written over max-plus semiring, 
where it is much easier to calculate. 

We remark that some aspects of the theory of attraction spaces have been investigated 
in [HI [iHl [30] in certain special cases. Also, the periodicity of max algebraic powers of 
matrices can be regarded from the viewpoint of max-plus semigroups as studied in |31j . 

2. Two TOPICS IN MAX ALGEBRA 

2.1. Spectral problem. Let A E M"^". Consider the problem of finding A G R+ and 
nonzero x G M.^ such that 



(1) 



A^ X = Xx. 
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If for some A there exists a nonzero x G M" which satisfies ([T]), then A is called a max- 
algebraic eigenvalue of A, and x is a max- algebraic eigenvector of A associated with A. 
With the zero vector adjoined, the set of max-algebraic eigenvectors associated with A 
forms a max cone, which is called the eigencone associated with A. 
The largest max-algebraic eigenvalue of A G Wl^"" is equal to 

n 

(2) A(A) = 0(TreA'=)^/^ 

k=l 

where Tr® is defined by Tif^^A) := 0"^^^ an for any A = {aij) G M"^". Further we explain 
the graph-theoretic meaning of ([2]), assumed that \{A) ^ 0. 

With A = (aij) G we can associate the weighted digraph V{A) = {N{A),E{A)), 

with the set of nodes N{A) = {1,. . . ,n} and the set of edges E{A) = {{i,j) \ a^j ^ 0} 
with weights w{i,j) = aij. Suppose that n = {ii, ip) is a path in V^A), then the weight 
of vr is defined to be w (vr. A) = ai^^i^ai^i^ . . . ai^_^i^ if p > 1, and 1 if p = 1. If ii = ip then 
TT is called a cycle. One can check that 

\{A) = max yu(cr. A), 
where the maximization is taken over all cycles in 'D{A) and 

denotes the geometric mean of the cycle cr = [ii, ...,ik,ii). Thus X{A) is the maximum 
cycle geometric mean of ^^(^4). 

A G M"^" is irreducible if for any nodes i and j there exists a path in 'P(A), which 
begins at i and ends at j. In this case A has a unique max-algebraic eigenvalue which 
equals A (A). 

Note that X{aA) = a\{A) and hence X{A/\{A)) = 1 if A(y4) > 0. Unless we need 
matrices with X{A) = 0, we can always assume without loss of generality that \{A) = 1. 
Such matrices will be called definite. 

An important relaxation of ([T]) is 

(3) A®x < Xx. 

The nonzero vectors x G M" which satisfy ([3]) are called subeigenvectors associated with 
A. With the zero vector adjoined, they form a max cone called subeigencone. This is a 
conventionally convex cone, meaning that it is closed under the ordinary addition. See 
[12] for more details. 

The eigencone (resp. subeigencone) of A associated with X{A) will be denoted by V{A) 
(resp. V*{A)). 
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2.2. Kleene stars. Let A E M"^". Consider the formal series 
(4) A* = I®A®A^®..., 

where / denotes the identity matrix with entries 

1, if ^ = i, 



0, otherwise. 



Series (jlj) is a max-algebraic analogue of {I — A) ^, and it converges to a matrix with 
finite entries if and only if X{A) < 1 [H [TT]. In this case 

(5) A* =I®A®A^®...® A''~\ 

which is called the Kleene star of A. 
For any A G 

(6) A is a Kleene star ■<=^ = A, an = 1 

The condition X{A) < 1 suggests that there is a strong interplay between Kleene stars 
and spectral problems. To describe this in more detail, we need the following notions and 
notation. 

A cycle a in is called critical, if fi{a,A) = X{A). Every node and edge that 

belongs to a critical cycle is called critical. The set of critical nodes is denoted by Nc{A), 
the set of critical edges is denoted by Ec{A). The critical digraph of A, further denoted 
by C{A) = {Nc{A) , Ec{A)) , is the digraph which consists of all critical nodes and critical 
edges of V{A). For definite A E it follows that aija*^ < 1 [3|. Further, 

(7) EE,{A)^aija*, = l. 

For definite A E M"^"", the relation between Kleene star, critical graph and spectral 
problems is briefiy as follows [31 [151 132] : 

(8) V*{A) = span(A*) = J a,A*„ a, E 

[ i=i 

(9) ViA) = I , a, E M 

[i€N,{A) 

(10) X E V*{A), {i,j) E Ec{A) aijXj = Xi. 

Equation ([8]) means that V*{A) is the max-algebraic column span of Kleene star A*, 
also called Kleene cone. This cone is convex in conventional sense. By ([9]), V{A) is the 
max subcone of V*{A), spanned by the columns with critical indices. Implication ( fTOi) 
means that for any subeigenvector x E V*{A) and i E Nc{A), the maximum in aijXj 
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is attained at j such that G Ec{A). In particular, {A x)i = Xi for all x G V*{A) 
and i G iVc(A). 

Not all columns in ([8]) and ([9]) are necessary. Let C{A) have ric G {1, . . . ,n} strongly 
connected components (s.c.c.) C^, for = 1, . . . , rzc. It follows from the definition of C{A) 
that s.c.c. are disjoint. The corresponding node sets will be denoted by A*"^. Let c 
denote the number of non-critical nodes of 'D{A). It can be shown [3l|15] that if z, j belong 
to the same s.c.c. of C{A), then the columns A* and A*- are multiples of each other. The 
same holds for the rows A* and A* . Hence 



where K is any set of indices which contains all non-critical indices and for every there 
is a unique index of this component in K. 

Consider A*j^j^, the principal submatrix of A* extracted from the rows and columns 
with indices in K. Condition ([6]) implies that A*j^j^ is itself a Kleene star. It follows 
from the maximality of that there is a unique permutation of K that has the greatest 
weight with respect to A*p^j^. The weight of a permutation tt of {1, . . . , n} with respect to 
A G W]^"^ is defined as nr=i '^ii^ii)- Thus A]^j^ is strongly regular in the sense of Butkovic 
[7j. From this it can be deduced that the columns of A* with indices in K are independent, 
meaning that none of them can be expressed as a max combination of the other columns. 
In other words [10], the columns of A* with indices in K (resp., in Nc{A) fl K) form a 
basis of V*{A) (resp., of V{A)). This basis is essentially unique [10], meaning that any 
other basis can be obtained from it by scalar multiplication. 

More precisely, the strong regularity of A\j^ is equivalent to saying that this basis is 
tropically independent, hence the tropical rank of A* is equal to ric + c, see [21 ETJ [28j| for 
definitions and further details. 



The transformation A ^ X ^AX is called a diagonal similarity scaling of A. Such 
transformations do not change \{A) and C{A) [21]. They commute with max-algebraic 



(12) 



(11) 




3. Visualization and Boolean matrices 



3.1. Visualization. Consider a positive x G and define 

/xi ... 0\ 



(13) X = diag(x):= : ••. ; 

VO ... xj 
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multiplication of matrices and hence with the operation of taking the Kleene star. Ge- 
ometrically, they correspond to automorphisms of R" , both in the case of max algebra 
and in the case of nonnegative linear algebra. Further we define scalings which lead to 
particularly convenient forms of matrices in max algebra. 
A definite matrix A G M"^" is called visualized, if 

(14) Uij < 1, Vi, j = 1, . . . ,n 

(15) a,, = 1, V(2,j)eE,(A) 

A visualized matrix A G M.^^^ is called strictly visualized if 

(16) aij = l^{t,j)eE,{A). 

Visualization scalings were known already to Afriat [1] and Fiedler-Ptak [23], and mo- 
tivated extensive study of matrix scalings in nonnegative linear algebra, see e.g. [2n [22| 
[361 [38]. We remark that some constructions and facts related to application of visualiza- 
tion scaling in max algebra have been observed in connection with max algebraic power 
method fl9\ [20] , behaviour of matrix powers |8] and max-balancing [36|, [38] . 

Visualization scalings are described in [12] in terms of the subeigencone V*{A) and its 
relative interior. For the convenience of the reader, we show their existence for any definite 
A G M"^"". In the proposition stated below, the summation in part 2. is conventional. 

Proposition 3.1. Let A G R"^" be definite and X = diag(x). 

1. If X = ^^^i A*^ then X~^AX is visualized. 

2. If X = Y17=i then X~^AX is strictly visualized. 

Proof. 1. Observe that x G V*{A) and x is positive. Then aijXj < Xi for all i,j implies 
xl^OijXj < 1, and by f[TU]) x^'^a^Xj = 1 for all (z, j) G Ec{A). 

2. Observe that x is positive, and that x G V*{A) since V*{A) is convex. Hence X~^AX 
is visualized. It remains to check that {i,j) ^ Ec{A) implies aijXj < Xt. We need to find 
k such that ajjO*^ < a*^. But this is true for /c = since a*j = 1 and aija*^ < 1 by ([7]). 
This completes the proof. □ 

More precisely [42], A G M"^" can be visualized by any positive vector in V*{A), and 
it can be strictly visualized by any vector in the relative interior of V*{A). 

3.2. Max algebra and Boolean matrices. Max algebra is related to the algebra of 
Boolean matrices. The latter algebra is defined over the Boolean semiring S which is the 
set {0, 1} equipped with logical operations "OR" a®h := a\J h and "AND" a®h := a /\h. 
Clearly, Boolean matrices can be treated as objects of max algebra, as a very special but 
crucial case. 
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For a strongly connected graph, its cyclicity is defined as the g.c.d. of the lengths of 
all cycles (or equivalently, all simple cycles) . If the cyclicity is 1 then the graph is called 
primitive, otherwise it is called imprimitive. We will not distinguish between cyclicity (or 
primitivity) of a Boolean matrix A and the associated digraph T>{A). Further we recall 
an important result of Boolean matrix theory. 

Proposition 3.2 (Brualdi and Ryser [6]). Let A G 5""^" be irreducible, and let be 
the cyclicity of V{A) (which is strongly connected). Then for each k > 1, there exists 
a permutation matrix P such that P~^A^P has r irreducible diagonal blocks, where r = 
gcd(/c,7^), and all elements outside these blocks are zero. The cyclicity of all these blocks 
is 'jA/r. 

In max algebra, let A e M"""". Define the Boolean matrix A^'^^ = (aj^^) by 



Let A,Be M"^". Assume that C{A) has Uc s.c.c. for n = 1, . . . ,nc, with cyclicities 7^. 
Denote by B^i, the block of B extracted from the rows with indices in A^^ and columns 
with indices in N,y. 

The following proposition can be seen as a corollary of Proposition I3.2[ The idea of the 
proof given below is due to Hans Schneider. See also [26] Section 3.1 and [8] Theorem 



Proposition 3.3. Let A e and X{A) ^ 0. 

1. A(yl^) = X'^iA). 

2. (A[^1)'= = (A'=)[^1. 

3. For each k > 1, there exists a permutation matrix P such that (P~^yl'^P)[^', for 
each /i = 1, . . . , ric, has := gcd(fc, 7^) irreducible blocks and all elements outside 
these blocks are zero. The cyclicity of all blocks in (P~^74'^P)|^' is equal to "^^/r^. 

Proof. We can assume that A is definite. Further, the diagonal similarity scaling com- 
mutes with max algebraic matrix multiplication and changes neither \{A) nor C{A) |21j . 
and by Proposition 13. part 2, there exists a strict visualization scaling. Hence we can 
assume that A is strictly visualized. In this case A^^^ = A'^l, where A'^l = {afh is defined 



(17) 




2.3. 



by 



(18) 
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It is easily seen that (At^l)*-' = {A'^^'^l As ^[^1 = A^'^\ all entries of A'^] outside the blocks 
aI^I are zero, which assures that (A^^^)^^ = {A\^}i)^. 

Proposition 13.21 implies that part 3. is true for (A'^^)'^ = (A'^Y^K This implies that 
P~^{A''y^^P has irreducible blocks and X{A'') = 1, which shows part 1. Also, P~^{A^y^^P 
has block structure where all diagonal blocks are irreducible and all off-diagonal blocks 
are zero. This implies (A'')^'^^ = (A'^Y^^ and parts 2. and 3. follow immediately. □ 

3.3. Cyclic classes. For a path P in a digraph G = {N,E), where N = {!,..., n}, 
denote by 1{P) the length of P, i.e., the number of edges traversed by P. 

Proposition 3.4 (Brualdi-Ryser [B]). Let G = {N,E) be a strongly connected digraph 
with cyclicity 7^. Then the lengths of any two paths connecting i E N to j E N (with i,j 
fixed) are congruent modulo 7^. 

Proposition 13.41 implies that the following equivalence relation can be defined: i ~ j if 
there exists a path P from i to j such that 1{P) = 0(mod7G). The equivalence classes of 
G with respect to this relation are called cyclic classes [H [391 SO]- The cyclic class of i 
will be denoted by [i]. 

Consider the following access relations between cyclic classes: [i] —^t [j] if there exists a 
path P from a node in [i] to a node in [j] such that 1{P) = t(mod7G). In this case, a path 
P with 1{P) = t(mod7G') exists between any node in [i] and any node in [j]. Further, by 
Proposition 13.41 the length of any path between a node in [i] and a node in [j] is congruent 
to t, so the relation [i] -^t [j] is well-defined. Classes [i] and [j] will be called adjacent if 

Cyclic classes can be computed in Od-El) time by Balcer-Veinott digraph condensation, 
where \E\ denotes the number of edges in G. At each step of this algorithm, we look for 
all edges which issue from a certain node and condense all end nodes of these edges 
into a single node. A precise description of this method can be found in [H |6]. We give 
an example of its work, see Figures 1 and 2. 




Figure 1. Balcer-Veinott algorithm 
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Figure 2. Result of the algorithm (left) and cyclic classes (right) 



In this example, see Figure 1 at the left, we start by condensing nodes 2 and 4, which 
are "next to" node 1, into the node 24. Further we proceed with condensing nodes 3 and 
5 into the node 35. In the end, see Figure 2 at the left, there are just two nodes 135 and 
246. They correspond to two cyclic classes {1,3,5} and {2,4,6} of the initial graph, see 
Figure 2 at the right. 

The notion of cyclic classes and access relations can be generalized to the case when G 
has ric disjoint components with cyclicities 7^, for /i = 1, . . . , ric (just like the critical 
graph in max algebra). In this case we write i ~ j if z, j belong to the same component 
and there exists a path P from i to j such that 1{P) = 0(mod7^). If 1{P) = t(mod7^), 
then we write [i] — [j]. In this case the cyclicity of G is 7 := 1cm 7^, /i = 1, . . . , 

We will be interested in the cyclic classes of critical graphs, and below we also give 
an explanation of these, in terms of the Boolean matrix A^'^'^. Let A G M"^". Following 
Brualdi and Ryser ^ we can find such ordering of the indices that any submatrix A 
which corresponds to an imprimitive component of C{A), will be of the form 



[C] 



(19) 



/ 






A 










[c\ 
S2S3 






A 






[c\ 

Sk-lSk 





where k is the number of cyclic classes in C^. Indices Sj and Sj+i for i = 1, . . . , — 1, 
and Sk and si correspond to adjacent cyclic classes. By Proposition 13.31 part 2, when A 
is raised to power fc, A^'"^ is also raised to the same power over the Boolean algebra. Any 
power of At"-^! has a similar block-permutation form. In particular, (A'^'')|15i' looks like 



(20) 











{A 





7m 



AC] 



••• 
••• 







\ 
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Theorem 5.4.11 of [29] implies that the sequence becomes periodic 

after k < {n — lY + 1, with period 7 = lcm(7^), /i = 1, . . . ,nc. In the periodic regime, 
all entries of nonzero blocks are equal to 1. 

4. Periodicity and complexity 

4.1. Spectral projector and periodicity in max algebra. We further assume that 
the critical graph C{A) occupies the first c nodes, i.e., that Nc{A) = {1, . . . , c}. 
For a definite A G M"^", consider the matrix Q{A) with entries 

c 

(21) = a*i^alj, i,j = l,...,n. 

k=l 

The max-linear operator whose matrix is Q{A), is a max-linear spectral projector associ- 
ated with A, in the sense that it projects M" on the max-algebraic eigencone {x \ A^x = 
x} [3]. We also note that Q{A) is important for the policy iteration algorithm of [13] . 
We will need the following property of Q{A) which follows directly from (I2T1) . 

Proposition 4.1. For a definite A G M"^", any column (or row) of Q{A) with index in 
1, . . . ,c is equal to the corresponding column (or row) of A* . 

This operator is closely related to the periodicity questions. 

Theorem 4.2 (BacceUi et al. [3J, Theorem 3.109). Let A G M"""" be irreducible and 
definite, and let all s.c.c. of C{A) have cyclicity 1. Then there is an integer T [A) such 
that A'' = Q{A) for all r > T{A). 

It can be easily shown that Theorem 14.21 can be generalized to the situation when A is 
in a blockdiagonal form with irreducible definite diagonal blocks An, . . . , A^^ so that 

Mil • • • \ 

(22) A= : ••. : 

V ... AuuJ 

If A is irreducible and definite, then A'' for all A; > 1 is in a blockdiagonal form (12^ . 
where all blocks An, . . . , A^u are irreducible and definite. 

Indeed, it is evident (when A is assumed to be visualised), that the m.c.g.m. of each 
block does not exceed 1. On the other hand, any column A* with 1 < z < c is a positive 
max-algebraic eigenvector of A and hence of A'^. This shows that the m.c.g.m. of each 
submatrix, being equal to the largest max-algebraic eigenvalue of that submatrix, cannot 
be less than 1. 

It follows from Proposition 13. 31 part 3 that all components of C(A'*') are primitive, where 
7 is the cyclicity of C{A). 
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These arguments lead us to the following extension of Theorem H? 

Theorem 4.3. Let A G M"^" be irreducible and definite, and let 7 be the cyclicity of 
C{A). There exists T{A) such that 

1. = Q(A^) for all t7 > T{A); 

2. A''+^ = A^ for all r > T{A). 

T{A) will be called the transient, and the powers A^ for r > T{A) will be called periodic 
powers. 

(r) 

It is also important that the entries alj , where i or j are critical, become periodic much 
faster than the non-critical part of A. The following proposition is a known result, which 
is proved here for convenience of the reader. 

Proposition 4.4 (Nachtigall [M])- Let A G M"^" be a definite irreducible matrix. Critical 
rows and columns of become periodic for r > n"^. 

Proof. We prove the claim for rows, and for columns everything is analogous. Let i G 
{1, . . . , c}. Then there is a critical cycle of length / to which i belongs. Hence a^^f'^ = 1 
for k > 1. Since for all m < k and any t = 1, . . . , n we have 

(ml) _ {{k-m)l) (ml) , (kl) 

it follows that 

(23) a(f)^0,M). 

m=l 

Entries a^-^'^ are maximal weights of paths of length k with respect to the matrix AK 
Since the weights of all cycles are less than or equal to 1 and all paths of length n are not 
simple, the maximum is achieved at k < n. Using (!23|) we obtain that a^f^^^''^ = afp for 
all t > n. Further, 

Hi+d) _ ri\(ti)(d) 



^is ' — ^ik ' ^ks ' 



and it follows that a'lf^^^'''^^^ = af^^^^^ for all t > ri and < d < l — l. Hence a^^^ is periodic 
for k > nl, and all these sequences, for any i = 1, . . . , c and any s, become periodic for 
k>n^. □ 



The number after which the critical rows and columns of A* become periodic will be 
denoted by Tc{A). 
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4.2. The ultimate spans of matrices. Max algebraic powers in the periodic regime 
have the following properties. 

Proposition 4.5. Let A G R"^"^ be a definite and irreducible matrix, and let t > be 
such that t'j > T{A). Then for every integer I > 

(24) A^^ = air^ ^r', A^r' = © «r , l<k<n. 

1=1 i=l 



Proof. Due to Theorem 14.31 for B = A^ and any r > T(B) we have 

c 

(25) 6M = 06*^6*^, l<k,j<n. 



i=l 



By Theorem 14.31 and Proposition 14.11 we have hi- = b^^^ = a^j^^^ and b*j = b^[^ = af^^ for 
all r > T{B) or equivalently t7 > T{A)., and any i < c. Hence 



(26) 47) = 047)a;f), l<k,j<n. 



,(i7) _ Jtj)Jtj) 

i=l 

In the matrix notation, this is equivalent to: 



(27) Al-, = ©alr^^r, A^2 = ©«lr^^*7, l<k<n. 

i=l 1=1 

Multiplying (p71) by any power A', we obtain (12^ . □ 
In the proof of the next proposition we will use the following simple principle 

(28) aJ)4J<a;r\ V2,j,/c,r,s, 
which holds for the matrix powers in max algebra. 

Proposition 4.6. Let A G M"^"" be a definite and irreducible matrix, and let i,j G 
{1, . . . , c} be such that [i] -^i [j], for some < Z < 7. 

1. For any r > Tc{A), there exists ti > such that 

(29) a^'^-'^Ar = ^5+^ a^'-^'^Al = 

2. If A is visualized, then for all r > Tc{A) 

(30) Ar = Al = 

Proof. If [i] —^i [j] then [j] — [i] where / + s = 7. By the definition of access relations 
there exists a critical path of length tij + / connecting i to j, and a critical path of 
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length + s connecting j to i. Hence afj^'^^'''^ a^^^"''^'^^ = 1, and in the visualized case 
aj*.i^+') = a5*/^+'^ = 1. Combining this with ([28]) we obtain 

J- 3' ^3 3'^ ^' J' 

Since r > n^, by Proposition KM Ar = ^^+(*i+*2+i)7 ^^^^ _ ^r+(ti+t2+i)7^ j^^^^^ 

inequahties (13T|) are equalities. Multiplying them by ctj-*^'^^'^ we obtain (!29l) . which is (l30ll 
in the visualized case. □ 

Proposition 14.61 says that in any power for r > n^, the critical columns (or rows) 
can be obtained from the critical columns (or rows) of the spectral projector QiA^) via 
a permutation whose cycles are determined by the cyclic classes of C{A). Proposition 14.51 
adds to this that all non-critical columns (or rows) of any periodic power are in the max 
cone spanned by the critical columns (or rows). From this we conclude the following. 

Proposition 4.7. All powers for r > T{A) have the same column span, which is the 
eigencone V{A'^) . 

Proposition 14 . 71 enables us to say that VlA"^) is the ultimate column span of A. Similarly, 
we have the ultimate row span which is ^((^1^)'^). These cones are generated by critical 
columns (or rows) of the Kleene star (A^)*. For a basis of this cone, we can take any set 
of columns {A'^)* (equivalently Q{A'^) or A"^ for r > T{A)), whose indices form a minimal 
set of representatives of all cyclic classes of C{A). This basis is tropically independent in 
the sense of [21 ESI EZl ■ 

4.3. Solving periodicity problems by square multiplication. Let A G M"^" and 
X{A) = 1. The t-attraction cone Attr(A, t) is the max cone which consists of all vectors 
X, for which there exists an integer r such that A^ ® x = A^^'^ ® x, and hence this is also 
true for all integers greater than or equal to r. Actually we may speak of any r > T{A), 
due to the following observation. 

Proposition 4.8. Let A be irreducible and definite. The systems A^ ^ x = ^4''+* x are 
equivalent for all r > T{A). 

Proof. Let x satisfy A^^x = A'*^*®^ for some s > T{A), then it also satisfies this system 
for all greater s. Due to the periodicity, for all k from T{A) < k < s there exists / > s 
such that A'' = AK Hence A'' (g) x = A''+^ (g) x also hold ior T{A) < k < s. □ 

Corollary 4.9. AttT{A,t) = Attr(A*, 1). 

Proof. By Proposition 14.81 Attr(A, t) is solution set to the system A^ ^ x = A''^* ® x for 
any r > T{A) which is a multiple of t, which proves the statement. □ 
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An equation of ® x = ^4''+* ® x with index in {1, . . . , c} will be called critical, and 
the subsystem of the first c equations will be called the critical subsystem. 

Proposition 4.10. Let A be irreducible and definite and let r > T{A). Then A"^ ® x = 
^ X is equivalent to its critical subsystem. 

Proof. Consider a non-critical component A^. ®a; = A^^*®a;. Using fl24p it can be written 
as 

(32) 04)A[.®a; = 04?A[.+*®x, 

i=l i=l 

hence it is a max combination of equations in the critical subsystem. □ 

Next we give a bound on the computational complexity of deciding whether x G 
Attr(y4, t), as well as other related problems which we formulate below. 

PI. For a given x, decide whether x G Attr(A,t). 

P2. For a given k : < /c < 7, compute periodic power A^ where r = A;(mod 7). 

P3. For a given x compute the ultimate period of {A^ x, r > 0}, meaning the least 

integer a such that A^'^" ® x = A*" ® x for all r > T{A). 

The following theorem is analogous to the results of Semancikova [39l |40] . 

Theorem 4.11. For any irreducible matrix A G R"^", the problems P1-P3 can be solved 
in 0{rt'\ogn) time. 

Proof. First note that we can compute both \{A) and a subeigenvector, and identify 
all critical nodes in no more than O(ra^) operations, which is done essentially by Karp 
and Floyd- Warshall algorithms |35]. Further we can identify all cyclic classes of C{A) by 
Balcer-Veinott condensation in 0{n'^) operations. 

By Proposition 14.41 the critical rows and columns become periodic for r > n'^. To know 
the critical rows and columns of a given power r' > T{A), it suffices to compute A"^ for 
arbitrary r > n"^ which can be done in O(logn) matrix squaring {A, A"^, A^, ...) and takes 
0(r;,^ log ra) time, and to apply the corresponding permutation on cyclic classes which 
takes O(n^) overrides. By Proposition 14. lOl we readily solve PI by the verification of the 
critical subsystem of A'"' ® x = A''^'^ ® x which takes 0{v?^ operations. Using linear 
dependence fl2^ the remaining non-critical submatrix of A!' , for any r > Ti^A) such that 
r = fc(mod 7), can be computed in Oiv?^ time. This solves P2. 

As the non-critical rows of A are generated by the critical rows, the ultimate period of 
^A" ® x\ is determined by the critical components. For visualized matrix we know that 
= A^'j. for all 1,3 such that [i] [j]- This implies = {A'®x)^ for [i] [j], 

meaning that, to determine the period we need only the critical subvector of A' ®x for any 
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fixed r > n^. Indeed, for any i G Nc{A) and r > ri^ tlie sequence ® t > 0} can 

be represented as a sequence of critical coordinates oi ®x determined by a permutation 
on 7^ cyclic classes of the s.c.c. to which i belongs. To compute the period, we take a 
sample of 7^ numbers appearing consecutively in the sequence, and check all possible 
periods, which takes no more than 7^ operations. The period of A^' ® x appears as the 
l.c.m. of these periods. It remains to note that all operations above do not require more 
than 0{n^) time. This solves P3. □ 

5. Properties of periodic powers 

5.1. Core matrix. In the sequel we always assume that A e R"^" is irreducible. Let 
C{A) consist of ric s.c.c. with cyclicities 7^, for // = 1, . . . , ric- Let 7 = l.c.m. (7^) and 
c be the number of non-critical nodes. Further it will be convenient (though artificial) 
to consider, together with these components, also "non-critical components" for ji = 
nc + 1, ■ ■ ■ ,nc + c, whose node sets consist of just one non-critical node, and the set of 
edges is empty. 

Consider the block decomposition of A^ for r > 1, induced by the subsets iV^ for 
H = 1, . . . ,nc + c. The submatrix of A"^ extracted form the rows in A*"^ and columns in N^, 
will be denoted by A^^^. If A is visualized and definite, we define the corresponding core 
matrix = (ct^,^), ix^u = 1, + c by 

(33) afj,^ = max{ajj- | i e iV^, j G A^,,}. 

The entries of (^A^°'^'^)* will be denoted by a*^. Their role is investigated in the next 
theorem. 

Theorem 5.1. Let A e R"^" he a definite visualized matrix and r > Tc{A). Let iJ>,Ty — 
1, . . . jTic + c be such that at least one of these indices is critical. Then the maximal entry 
of the block A^^^ is equal to a*^^. 

Proof. The entry q;*^, is the maximal weight over paths from n to z/, with respect to the 
matrix A"-^"^^. We take such a path (/xi, . . . , ni) with maximal weight, where /ii := fi and 
Hi ~ V. With this path we can associate a path tt in V[A) defined by tt = ti o ai o r2 o 
. . . o ai-i oTi, where are critical paths which entirely belong to the components C^., and 
(7j are edges with maximal weight connecting C^. to C^-^-^. Such a path exists since any 
two nodes in the same component can be connected to each other by critical paths 
if II is critical, and if /i is non-critical then consists just of one node. The weights of 
Ti are equal to 1, hence the weight of tt is equal to q;*^. It follows from the definition of 
a^i, and a*^^ that q;*^ is the greatest weight over all paths which connect nodes in to 
nodes in Ci,. As at least one of the indices fi, v is critical, there is freedom in the choice 
of the paths t\ or ti which can be of arbitrary length. Assume w.l.o.g. that is critical. 
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Then for any r exceeding the length of ai o r2 o . . . o cr^^i o n which we denote by /^j^, the 

(r) 

block Altu contains an entry equal to a*^^ which is the greatest entry of the block. Taking 
the maximum T'{A) of l^j^^, over all ordered pairs (/x, z/) with fi or u critical, we obtain the 
claim for r > T'{A). Evidently, T'{A) can be replaced by Tc{A). □ 

5.2. CSR-representation. For a definite visualized matrix A G M"^", the statements of 
Propositions 113] and HJ] can be combined in the following. Let C E M"^'' and R G 1^+^" 
be matrices extracted from the first c columns (resp. rows) of Q{A'^) (or equivalently 
{A"')*), and let S := A^'^\ the critical matrix of A defined by (ITTl). 

Theorem 5.2 (Schneider [37]). Let A G M"^" be definite and visualized. For r > T{A) 
and r = /(mod 7), = C O S" i?. 

Proof. By (Ell) Q{A^) =C®R, and by Theorem HJl^^* = Q{A^) = for 7^ > T{A). 

Thus the claim is already proved for r = 'yt > T{A). 

Note that C (resp. R) can be extracted from the first c columns (resp. rows) of A'^^. 

As S = A^'^\ it follows that 1) the {i,j) entry of 5' can be 1 only if [i] [j], 2) for 
each pair of classes [i] [j] and each ii G [i] there exists ji G [j] such that the («i, ji) 
entry of 5' equals 1. 

Using these two observations and equation (129!) applied to A^'^^^ and for jt > 

T(A), we obtain that the critical columns of A^'^^'^''^ are given by C (S> 5*' and the critical 
rows of y4 are given by S'' ® R. 

Combining this with any of the two equations of (El]), we obtain that A'^*"''' = C®S''®R, 
for any 7* > T(v4). □ 

Further we observe that the dimensions of periodic powers and the CSR-representation 
established in Theorem 15.21 can be reduced. 

The rows and columns with indices in the same cyclic class coincide in any power 
y4'', where r > Tc{A) and A is definite and visualized. Hence the blocks of , for 
/i, z/ = 1, . . . , ric + c and r > Tc{A)., are of the form 

(34) 4;) 



where k (resp. m) are cyclicities of (resp. Cj,), indices (resp. t\,...tm) 

correspond to properly arranged cyclic classes of (resp. Cy)., and E^.i. are matrices 
with appropriate dimensions with all entries equal to 1. We assume that has just one 
"cyclic class" if /i is non-critical. 
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defines the matrix i^^') G >^ (^'+^) 



, where c is the total number of cyclic 



classes, as matrix with entries a^g"].. Corresponding to (IMl) . this matrix has blocks 



(35) 



"'siti 



r) \ 



~{r) 



J 



It follows that = v4(''i)(g)v4(''2) for all ri, > Tc(A). In words, the multiphcation 

of any two powers A^^^'^^ and A^^'^^ for ri, r2 > Tc(y4) reduces to the multiplication of A^'^^'^ 
and 

If we take r = 'jt + I > T(A) (instead of Tc{A) above) and denote A := A^'^^'^^\ then 
due to the periodicity we obtain 



(36) 



so that A^^'^ can be regarded as the rth power of A, for all r > T{A). 

Matrices C, R, and S** for t > (n — 1)^ + 1 have the same block structure as in 
This shows that the behavior of periodic powers A^ is fully described by 

(37) = CS^R, r = /(mod 7), 

where is a c x c Boolean matrix with blocks 



(38) 



MM 



/o 1 

1 



\1 



o\ 



1 

0/ 



Sf,^ = for /i 7^ I/, 



and C e m(^^+^)^^ and R G mJ'(^+^) are extracted from critical columns and rows of 
A'^ = Q{A'^), or equivalently, formed from the scalars in the blocks of C and R. 

5.3. Circulant properties. Matrix A G M"^" is called a circulant if there exist scalars 
a^i such that Qij = ad whenever j — i = d{mod n). This looks like 

( a^ a'} a% ■ ■ ■ an \ 
(39) A = an-i an ai ■ ■ ■ an-2 



We also consider the following generalizations of this notion. 
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Matrix A G M™^" will be called a rectangular circulant, if aij = aps whenever p = 
i + t(mod m) and s = j + t(mod n), for all i,j,t. When m = n, this is an ordinary 
circulant given by fl39|) . 

Matrix A G M™^" will be called a block kxk circulant when there exist scalars ai, . . . ,ak 
and a block decomposition A = {Aij), i,j = 1, . . . ,k such that Aj = ctdEij if j — i = 
d{mod k), where all entries of blocks Eij are equal to 1. 

A G M™^" is called d-periodic when = if (s — j) modn is a multiple of d, and 
when Qji = Qsi if {s — j) modm is a multiple of rf. 

We give an example of 6 x 9 rectangular circulant A: 

/O 1 2 1 2 1 2\ 
2 1 2 1 2 1 
^_ 120120120 
012012012 
2 1 2 1 2 1 
^1 20120120y 

This example provides evidence that a rectangular mxn circulant consists of ordinary 
d X d circulant blocks where d =g.c.d.(m, n). In particular, it is d-periodic. Also, there 
exist permutation matrices P and Q such that PAQ is a block circulant: 



B = 



(o 








1 


1 


1 


2 


2 


2\ 











1 


1 


1 


2 


2 


2 


2 


2 


2 











1 


1 


1 


2 


2 


2 











1 


1 


1 


1 


1 


1 


2 


2 


2 











V 


1 


1 


2 


2 


2 








0/ 



We formalize these observations in the following. 

Proposition 5.3. Let A G M'J^^" be a rectangular circulant, and let d =g.c.d.{m,n) . 

1. A is d-periodic. 

2. There exist permutation matrices P and Q such that PAQ is a block dxd circulant. 

Proof. 1. There are integers ti and t2 such that d = tim + t2n. Using the definition 
of rectangular circulant we obtain a^j = a^s, if s = j + tim(mod n), and hence if s = 
j + d{mod n). Analogously for rows, we obtain that aji = asi, if s = j + t2n{mod m), 
and hence if s = j + d{mod m). 

2. As A is (i-periodic, all rows such that i+d = j(mod m) are equal, so that {1, . . . , m} can 
be divided in d groups each with m/d indices, in such a way that Ai. = Aj. if i and j belong 
to the same group. We can find a permutation matrix P such that A' = PA will have 
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rows A[, = . . . = A'^, = Ai., A'^^^, = . . . = A'2^, = A2., and so on. Analogously we can find 
a permutation matrix Q such that A" = PAQ will have columns A',[ = . . . = A"^ = A'.^, 
A"^_^_^ = . . . = = A',2, and so on. Then A" has blocks (A'-j) for i,j = l,...,d. of 
dimension n/d x m/d, where A'-j = aijEij, and Eij has all entries 1. As A is li-periodic, 
it can be shown that its submatrix extracted from the first d rows and columns, is a 
circulant. Hence A" is a block circulant. □ 

When g.c.d.(m, n) = 1, the rectangular circulant is 1-periodic, i.e., a constant matrix. 

Proposition 5.4. Let A G M"^" be a definite visualized matrix which admits block de- 
composition ( |34|) . and r > T{A). Let C^,Cy be two (possibly equal) components of C{A), 
and d = g.c.d.{;~^^, •y^). 

1. A^jI} is a rectangular circulant (which is a circulant if fi = u). 

2. For any critical fi and v, there is a permutation P such that {P^AP)^^i^ is a block 
d X d circulant matrix. 

3. If r is a multiple 0/7, then A^^l are circulant Kleene stars, where all off-diagonal 
entries are strictly less than 1. 

Proof. 1.: Using Eqn. (l30l) we see that for all [i,]) and {k, I) such that k = i + t(mod7^) 
and I = j + t(mod7i/), 

Sfct; ^itl ^itj 

2. : If n = u then P = I, and if fi ^ u then P is any permutation matrix such that its 
"subpermutations" for A*"^ and N^, are given by P and Q of Proposition 15.31 

3. : Part 1 shows that A^^^ are circulants for any r > T{A) and critical /i. If r is a multiple 
of 7, then A^^l are submatrices of = Q{A'^) and hence of {A'^)*. This implies, using 
(E]), that they are Kleene stars. As the fith component of C{A) is just a cycle of length 7^, 
the corresponding component of C{A'^) consists of 7^ loops, showing that the off-diagonal 
entries of Afj,ji are strictly less than 1. □ 

6. Attraction cones 

6.1. A system for attraction cone. Let A E R"^" and \{A) = 1. Recall that attraction 
cone Attr(y4, t), where t > 1 is an integer, is the set which consists of all vectors x, for 
which there exists an integer r such that A^ ® x = A^^^ ® x, and hence this is also true 
for all integers greater than or equal to r. We have shown above, see Subsect. 4.3, that 
Attr(y4, t) is solution set to the critical subsystem of the system A^^^^x = A^ ^x, for any 
f > Tc{A). Next we show how the specific circulant structure of A"^ at r > Tc{A) can be 
exploited, to derive a more concise system of equations for the attraction cone Attr(A, 1). 
Due to Theorem 15.11 the core matrix Af^°'^'^ = {a^^ \ fi,^ = l,...,nc}, and its Kleene 
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Star {A^"''^y = {a*^ \ fi,u = 1, . . . ^n^ } will be of special importance. We introduce the 
notation 

(40) 

ir«(^) = {t>c I a,? = zeN„ 

where and Ci, are s.c.c. of C{A), and Ni, are their node sets, and uit) in the second 
definition denotes the index of the non-critical component which consists of the node t. 
The sets M^.^'^i) are non-empty for any r > Tc{A), due to Theorem 15. 1[ 

The results of Subsect. 15.31 lead to the following properties of Mu\i) and K^''\i). 

Proposition 6.1. Let r > Tc{A) and /i, G {1, . . . , n^}. 

1. // [z] [j] and I, J e then M^"+*^(i) = M!f\j) and K^''+'\i) = K^''^]). 

2. Each My^\i) is the union of some cyclic classes ofC^. 

3. Let i G Nfj^ and d = (7.0.^.(7^,71,). Then, if [p] C My '\i) and [p] —>-d [s] then 
[s] C M!)'\i). 

4. Let i,j e andp,s e Ny. Let [i] [j] and [p] [s\- Then [p] C Ml''\i) if 
and only if [s\ C My''\j). 

Next we establish the cancellation rules which will enable us to write out a concise 
system of equations for the attraction cone Attr(74, 1). 
If a < c, then 

(41) {x: ax®h = cx®d} = {x:b = cx®d}. 
Now consider a system of equations over max algebra: 

n n n 

(42) ^ aiiXi © ci = ^ a2iXi © C2 = . . . = ^ a„jXi © c„. 

i=l i=l i=l 

Suppose that ai, . . . , a„ G IR+ are such that an < at for all I and i, and Si = {i \ an = a,} 
for Z = l,...,n. Let Si be such that [j^^i Si = {l,...,n}. Repeatedly applying the 
elementary cancellation law described above, we obtain that (l42l) is equivalent to 

(43) ^ aiXi © ci = ^ aiXi © C2 = . . . = ^ aiXi © c„. 

We will refer to the equivalence between fH2l) and fH3|) . which we acknowledge to [18], as 
to chain cancellation. 

Using notation (140!) and Proposition 16.11 we formulate the main result of the paper. 
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Theorem 6.2. Let A G M"^" be a visualized matrix and r > T{A) be a multiple 0/7. 
Then the system A'^ ® x = A^~^^ ® x is equivalent to 

© < J xA® al,^^t)Xt = 

(44) " 

= Xk © I Xk j © Oi*^i,{t)Xt-, 

where is the component of C{A) which contains both [i] and [j], and [i] and [j] range 
over all pairs of cyclic classes such that [i] ^1 [j]. 

Proof. By Proposition 14.101 A*" ® x = A"^^^ © x is equivalent to its critical subsystem. 
Consider a critical component of © x = A"^^^ © x: 

(45) 0aj)x, = 0ar')x„ . = l,...,c. 

k k 

Consider j such that [i] ^1 [j]. Then by Proposition 14.61 a^l^^^ = a^-)^, hence the critical 



"jk ' 

subsystem of A'' © x = A^'^^ © x is as follows: 

(46) 0a;^x, = 04Jx„ W^,J : [z] [j]. 

k k 

The number of non-identical equalities in fH5l) and fj46|) is equal to the total number of 
cyclic classes. 



Proposition 15.41 part 3, implies that all principal submatrices of A^ extracted from 
critical components have a circulant block structure. In this structure, all entries of the 
diagonal blocks are equal to 1, and the entries of all off-diagonal blocks are strictly less 
than 1. Hence we can apply the chain cancellation (equivalence between ( H2l) and ( l43i) ) 
and obtain the first terms on both sides of f j44l) . By Theorem 15.11 each block A^i, contains 
an entry equal to a*^. For a non-critical i^(t), this readily implies that the corresponding 
"subcolumn" A^y{t) contains an entry ^^^.(j)- Applying the chain cancellation we obtain 
the last terms on both sides of (jllj). Due to the block circulant structure of A^y with 
both n and v critical, see Proposition 15.41 or Proposition 16.11 we see that each column of 
such block also contains an entry equal to a*^. Applying the chain cancellation we obtain 
the remaining terms in ( 14^ . □ 

As Attr(y4,t) = Attr(74*, 1), system (jHj) also describes more general attraction cones, 
it only amounts to substitute C{A^) for C{A) and the entries of ((A*)'^'""^)* for a*^ (the 
dimension of this matrix will be different in general, see Proposition 13.31 part 3). 

We note that the system for Attr(A, 1) naturally breaks into several chains of equations 
corresponding to the s.c.c. of C{A). If we start with ( H6i) . it can be equivalently written as 
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R® X = H ®y where R is the factor in CSi?- representation, and H G M^^"" is a Boolean 
matrix with entries 



(47) /i, 



1, if^eiV^, 

0, otherwise. 



We can apply cancellation as described in the proof of Theorem l6.2l to get rid of redundant 
terms on the left-hand side of the two-sided system. 

If C{A) is strongly connected then H is a, vector of all ones, and the two-sided system 
R®x = H ®y becomes essentially one-sided. We treat this case in the next subsections. 

6.2. Extremals of attraction cones. System (jUj) in general consists of several chains 
of equations corresponding to s.c.c. of C{A). Each chain is of the form 

(48) OjXj = O-iXi = . . . = O^iXi, 

where Ti U . . . U = {1, . . . , n} and come from the entries of {A^°'^'^)*. 

Here we consider only the case of strongly connected C{A), i.e., only one chain. By 
scaling Ui = aiXi we obtain 

(49) ®y^ = ®y^ = ■■■ = ®y. 

Note that when the critical graph is not strongly connected, we have several chains 
of equations and the coefficients of psj) in general cannot be scaled to get fH9l) for each 
chain at the same time. 

By e* we denote the vector which has the ith coordinate and all the rest equal to 0. 
Vector y e will be called scaled if 0"=]^ i/i = 1, and set S is called scaled if it consists 
of scaled vectors. We say that V C is generated by S C (also, is a generating 
set of y) if 5 C and for any x & V there exist y^, . . . ,y'' E S and ai, G M+ such 
that X = aiy\ 

We investigate extremal solutions of fH9|) : a solution x is called extremal ii x = y®z for 
two other solutions y, z implies that x = y oi x = z [10]. The following can be deduced 
from the results of P 110] • 

Proposition 6.3. The solution set of any finite system of max-linear equations has a 
finite generating set. In particular, it is generated by extremal solutions, and any set of 
scaled generators for the solution set contains all scaled extremal solutions. 



In the next proposition we show that extremal solutions of (H9!) can have only and 1 
components. 
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Proposition 6.4. Let y be a scaled solution of and let < yi < 1 for some i. Then 
y is not an extremal. 

Proof. Let := {i \ < y^ < 1} and := {i \ yi = 1}, and define vectors and v^{k) 
for eacli k G by 

(50) 'f'^'^", 'f'^*-'^W, 

I 0, otlierwise I 0, otlierwise 

Observe tliat botli f° and v^{k) for any k G -ft'"^, are solutions to fH9l) . different from y. 
We fiave: 

(51) y = t;0© 

fience ?/ is not an extremal. □ 
Let T = [tij) be the m x n — 1 matrix defined by 

0, otherwise, 



(52) t 



where Tj are from fl49l) . 

A subset K C {1, . . . , n} is called a covering of T if each Tj contains an index from K. 
The following is immediate. 

Proposition 6.5. A scaled vector y is a solution of (HOj) if and only if := {i\yi = 1} 
is a covering ofT. 

A covering K is called minimal if it does not contain any proper subset which is also a 
covering. 

A covering K will be called nearly minimal if it contains no more than one proper 
subcovering K' . Observe that then the complement K\K' consists of just one index. 
Hence, a covering is nearly minimal if and only if there may exist no more than one i & K 
such that K\{i} is a covering. 

Proposition 6.6. Extremal solutions of fHOl) are precisely = ^^g/^ej, where K is a 
nearly minimal covering of T. 

Proof. If a covering K is not nearly minimal, then there exist i and j such that K[i] : = 
K\{i} and K[j] := K\{j} are both coverings of T. Then f ■'^W y^U] are both solutions 
and V = f © ^^[i] hence is not extremal. 

Conversely, if is not extremal, then there exist y ^ and z ^ such that 
= y ® z. Evidently y < and z < . By Proposition 16.41 we can represent y and z 
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as combinations of — 1 solutions of ( l49l) . These solutions correspond to coverings, which 
must be proper subsets of K. At least two of these coverings must be different from each 
other, hence K is not nearly minimal. □ 

Thus, the problem of finding all nearly minimal coverings of a — 1 matrix is equivalent 
to the problem of finding all extremal solutions of (149|) . 

The following case applies if the critical graph is strongly connected and occupies all 
nodes. 

Corollary 6.7. Let Ti, . . . , he pairwise disjoint, then the scaled extremals are precisely 
all vectors = ©jg^e*, where S is an index set which contains exactly one index from 
each set Ti. 

Proof. Any such set S forms a minimal covering of T, and it can be shown that the 
solution set of (H9|) is generated by , so there are no more scaled extremals (or nearly 
minimal coverings). □ 



6.3. An algorithm for finding the coefficients of an attraction system. Coeffi- 
cients of the system of equations which defines attraction cone are determined by the 
entries of (^A'""'^'^)* which can be found in 0{{nc + c)^) operations, where Uc is the number 
of s.c.c. of C{A) and c is the number of non-critical nodes. However it remains to find the 
places where these coefficients appear, i.e., the sets M!.''\i) and K^'^^i) for z = 1, . . . , c. 
Solving this problem, we get another polynomial method for computing the coefficients 

of dSD. 

Here we restrict our attention to the case when C{A) is strongly connected. In this case 
there are no second terms on both sides of (jS]) and we need only K^^\i). The digraph 
2)(^^Core-j associated with the matrix A^°^'^ consists of one critical node which corresponds 
to the whole C{A) and will be denoted by /x, and c non-critical nodes z/(t), for t > c. The 
core matrix is of the form 

(53) = 




the entries a^^, and the entries of h, g and B = {K(s),u{t)) are given by 
(54) /i^{t) = = max afct, g^it) = a^it)^ = maxatk, t > c, 

k=l k=l 
K{s)u{t) = (^i,(^s)u{t) = Cist, S > C, t > C. 
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Denote by [—*m i] the cychc class [j] such that [j] 
initiaUze Boolean c-vectors Pt by 



[i]. For each t 



1, . . . , c, we 



(55) 



PS) 



if 

otherwise. 



1 z] n arg max akt 

k=l 



Pt encode the Boolean information associated with h (shifted cyclically by 1). 

Further we compute the Kleene star of the non-critical submatrix B := Amm, where 
M denotes the set of non-critical nodes, and store the information on the lengths of paths 
with maximal weight and length not exceeding c in sets Ugt associated to each entry of 
B. To compute these sets we use the formula 

(56) B* = I®B®...®B~'-\ 



where for each entry of this matrix series we find the arguments of maxima. 

To combine the information associated with h and B*, we recall the max- algebraic 
version of bordering method [H], which computes 



(57) 



C\* 




B* ®g B* 



h'^^B* 
B* ®g®h^ 



B* 



where h,g E ]R^j_. Note that all information that we need for system (jH]), is in the entries 
of h'^ (8> B* and in the indices of equations of the system where the entries of h'^ ® B* 
appear. Computing (/i^ ® B*)i means in particular obtaining the "winning" indices 



(58) 



Wt = argmax/ii.(s)6* 



s>c 



After that, the idea is to combine Ps with Ust for all s E Wt and unite the obtained 
indices. More precisely, for each number m stored in Ust we define the shifted Boolean 
vector P^^ by 



(59) 

Equivalently, 
(60) 

After that, we define 
(61) 



if [- 
otherwise. 



1 i] n arg max aks 7^ 

k=l 



p: 



1 ^ Ps(j) = 1 and [j] [i 



Gr.= \/ V p: 

s£Wt m&Ust 



Proposition 6.8. Let t > c, and let r > T{A) be a multiple of Then for all i < c, 
t e ir(")(z) z/ and only if Gtii) = 1. 
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Proof. Gt{i) = 1 if and only if there exist s ^ Wt and m G Ugt such that P^''^{i) = 1. 
Then there exists a path tti o r o 7r2 where tti starts in i, belongs to C{A) and has length 
— m — l(mod 7), r is an edge which attains /ij,(<j) = max^^^a/s, and 7T2 is entirely non- 
critical, has length m, weight Vl(^s)v(t)^ ^^"^ connects s to t. The weight of tti o r o is 
equal to 1 ■ /i^ ■ h*^ = and the length is a multiple of 7, meaning that a-[^ = 

and t G 

The other way around, let t G K^'''\i). Then there exists a path P of length r that 
connects z to t and has weight = ^i^(s)K{s)u{t)- '^^^ decompose P = n'lOr'oni^, 

where r' is an edge connecting a critical node to a non-critical node, and vTg has only non- 
critical nodes. Obviously w{t' 01x2) < 0^ hy(s)bli^s)u{t) = ^IluH)- w{P) = hence 
w{7t[) = 1 so that tt[ entirely belongs to C{A), and w{t' o ti!^) = Q^^jy(t)- In particular, 
w(r') = h^(^s) and ^(vrg) = b*^(^^y^^-^ for certain s > c. From this we conclude that Gtii) = 
1. □ 

Summarizing above said, we have the following algorithm for computing the coefficients 
of (l44l) in the case when C{A) is strongly connected. Recall that in this case there is no 
second term on both sides of (jSj). The computation of coefficients of the third term 
includes the computation oi B* and the sets K^"^"* (i) for i < c (in fact we can improve 
the algorithm since only 7 of them are different). 



ALGORITHM 1. Compute the coefficients of f j44l) if C{A) is strongly connected. 

Input. Visualized matrix A, critical graph C{A) which is strongly connected and the 
cyclic classes of C{A). 

1. Compute h and initialize Pt for t > c. This takes cc operations. 

2. Compute B* and sets Ust for all s,t > c. It takes O(c^) operations (see (l56ll ). 

3. Compute hB* and Gt for t > c, by ([58]), (EQl) and (Ell). Computation of hB* and Wt 
by (158|) requires 0{(f) operations, computation of shifted Boolean vectors is O(cc^), and 
the conjunction (I6T!) takes 0{cc^) operations. 

4. Compute K'^^\i) using Proposition 16.81 This requires cc operations. 

As the overall complexity does not exceed 0{(f)+0{c(^) = 0{nc^) operations, we conclude 
the following. 

Proposition 6.9. Let A G M"^" be visualized, C{A) be strongly connected, c be the number 
of non-critical nodes, and suppose we know C{A) and all cyclic classes. Then Algorithm 
1 computes the coefficients of the attraction system in no more than 0{7?n) operations. 

It is also important that the eigenvalue and an eigenvector of irreducible matrix can be 
computed by the policy iteration algorithm of [12], which is very fast in practice. After 
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that, C{A) and the cychc classes can be computed in O(n^) time. Thus we are led to 
an efficient method of computing the coefficients in the case when A is irreducible and 
C{A) is strongly connected, especially in the case when the number of non-critical nodes 
is small. Note that the case of irreducible A and strongly connected C{A) is generic when 
matrices A are real and generated at random. Also, in this generic case it almost never 
happens that maxima in blocks or among the weights of paths are achieved twice, which 
means that we do not need to assign Boolean vectors or sets to each entry. In this case 
the total number of operations in the algorithm is reduced to O(c^) + 0{cc). 

When C{A) is not strongly connected, the bordering method ( l57j) can be used to obtain 
an algorithm which operates only with the entries of A'""^'^. However, the complexity of 
operations with indices and Boolean numbers significantly increases in that general case. 

7. Examples 

7.1. Square multiplication. In this subsection we will examine the problems that can 
be solved by matrix squaring on 9 x 9 real matrix over the max-plus semiring: 



/-I 





-1 


-1 


-9 


-7 


-10 


-4 


-8\ 





-1 





-1 


-10 


-1 


-10 


-9 


-4 


-1 


-1 


-1 





-2 


-3 


-2 


-6 


-6 





-1 


-1 


-1 


-10 


-6 


-10 


-6 


-1 


-10 


-2 


-8 


-1 


-1 





-1 


-10 


-1 


-5 


-5 


-10 


-9 


-1 


-1 





-3 


-6 


-9 


-10 


-7 


-10 





-1 


-1 


-8 


-8 


-75 


-80 


-77 


-83 


-80 


-77 


-82 


-2 


-0.5 


^-84 


-81 


-77 


-80 


-78 


-77 


-78 


-0.5 


-2 J 



The corresponding max-times example is obtained by, e.g., taking exponents of the 
entries. 

The critical graph of A, see Figure 3, has two s.c.c: Ci with nodes Ni = {1,2,3,4} 
and C2 with nodes N2 = {5, 6, 7}. The cyclicity of Ci is 71 = 2 and the cyclicity of C2 is 
72 = 3, so the cyclicity of C{A) is 7 = lcm(2, 3) = 2 x 3 = 6. 



6 




1 2 ^ 3 ^ 4 5 < 7 



Figure 3. The critical graph of A 
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The matrix can be decomposed into blocks 



A 



/ Au Ai2 Aim \ 

A21 A22 A2M 
\Ami Ami Amm / 



where the submatrices An and A22 correspond to two s.c.c. Ci and C2 of C{A), see Figure 
3. They equal 



An 



/-I -1 -l\ 
0-10-1 

-1 -1 -1 

\0 -1-1 -1/ 

and Amm is the non-critical principal submatrix 



A^ 



22 




A 



MM 



-2 -0.5 
-0.5 -2 



The submatrices A12, A21, Aim and A2M are composed of randomly taken numbers from 
— 1 to —10, and Ami and Am2 are composed of randomly taken numbers from —75 to 
-85. 

It can be checked that the powers of A become periodic after T{A) = 154. 

We will consider the following instances of problems P2 and P3. 
P2. Compute A"^ for r > T{A) and r = 2(mod 6). 
P3. For given x G M^, find ultimate orbit period of A'^ ® x. 

Solving P2. Using the idea of Proposition I4.11[ we perform 7 squarings A, A"^, A"^, . . . 
to raise A to the power 128 > 9 x 9. This brings us to the matrix 



/.(128) .(128) 4(128)\ 



A 



128 



HI 

1(128) 
^21 



^2 

,(128) 
^22 



IM 

,(128) 
^2M 



. .(128) .(128) .(128) 
\^M1 ^M2 ^MM , 



where 



4(128) 
^11 



all entries of A 



(128) 
12 



/o 

-1 



V-1 

,(128) 



-1 -l\ 
0-10 
-1 -1 
-1 0/ 



A 



(128) 
22 







1 o\ 

1 -1 

-V 



and ^21^°^ are —1 and 
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A 



(128) 
IM 



A 



(128) 
Ml 



v- 



/- 


-2 


5 -1\ 




-1 


5 -2 






-2 


5 -1 




V- 


-1 


5 -2/ 


76 




-75.5\ 


T 


75 




-76.5 




76 




-75.5 




75 




-IQ.hj 





A, 



(128) 
2M 



(128) 
M2 



/-1.5 
-2.5 
\-2.5 



-2\ 
-2 



76 -76.5\ 
76 -76.5 
y-76 -76.5y 



We are lucky since 128 = 2 (mod 6), as we already have true critical columns and rows 
of . However, the non-critical principal submatrix of A^'^^ is 



A 



(128) 
MM 



-64 -65.5 
-65.5 -64 



It can be checked that this is not the non-critical submatrix of A"' that we seek (recall that 
T{A) = 154). Hence, it remains to compute the principal non-critical submatrix A^^j^. 

We note that A^^'^ has critical rows and columns of the spectral projector Q{A), since 
132 is a multiple of 7 = 6. In A^^"^, the critical rows and columns 1 — 4 (in Ci) are the 
same as that of A^^^, since 71 = 2 and both 128 and 132 are even. The critical rows 5 — 7 
(in C2) can be computed from those of A^"^^ by cyclic permutation (5, 6, 7) — > (7, 5, 6), and 
the critical rows 5 — 7 can be computed by the inverse permutation (5, 6, 7) — » (6, 7, 5). 
This implies that all blocks in A^"^^ are the same as in A^'^^ above (in the analogous block 
decomposition of A^^'^), except for 



A. 



(132) 
22 



/o 

-1 

V-1 



4 



-i\ 
-1 

0/ 



A. 



(132) 
2M 



/-2.5 
-2.5 



-2\ 
-1 

-y 



Now the remaining non-critical submatrix of A"^ can be computed using linear depen- 
dence (^^, which specifies to 



1=1 



This yields 



A 



(r) 
MM 



-76.5 -77 
-78 -76.5 
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Solving P3 We examine tlie orbit period of A^x for X X J Ob ^ Ob J Ob , where 







= [1 


2 3 


4 


5 


6 


7 


8 


9] 








x"^ 


= [1 


2 3 


4 














0] 








x^ 


= [0 


1 


1 








1 


1 


1] 








4 

X 


= [0 


1 


1 














0] 








for X = 


1 2 'H 


4 

X . 




















(g) X^ 


= [8 


7 


8 


7 


7 


7 


8 


X 


X 




^^128 


x"^ 


= [3 


4 


3 


4 


3 


3 


3 


X 


X 


y3 


^^128 


(g) X^ 


= [1 


1 


1 


1 


1 








X 


X 




^^128 


® x^ 


= [1 


1 


1 


1 











X 


X 



Here x correspond to non-critical entries which we do not need. The cyclic classes of Ci are 
{1,3}, {2,4}, and the cyclic classes of C2 are {5}, {6} and {7}. From the considerations 
of Proposition 14.111 it follows that the coordinate sequences {{A'^x)i, r > T{A)} are 

2/1, 2/2, yi, 2/2, • • • , for i = 1,2,3,4, 

2/5, 2/6, 2/7, 2/5, 2/6, 2/7, • • • , for i = 5, 6, 7. 

Looking at y^,. . . above, we conclude that the orbit of x^ is of the largest possible 
period 6, the orbit of x'^ is of the period 2 (in other words, G Attr(A, 2)), the orbit of 
x^ is of the period 3 (i.e., x^ G Attr(A, 3)), and the orbit of x^ is of the period 1 (i.e., 
x^ G Attr(A, 1)). 

7.2. Circulants. Here we consider a 9 x 9 example of definite and visualized matrix in 
max-plus algebra 



(62) A = 



(-8 





-1 


-8 


-8 


-9 


-4 


-5 


-l\ 


-4 


-5 





-2 


-6 





-7 


-3 


-9 


-7 


-9 


-8 





-8 


-4 


-6 


-9 


-10 


-8 


-8 


-10 


-7 





-4 


-6 


-10 


-1 


-2 


-8 


-7 


-4 


-8 





-3 


-1 


-10 





-1 


-2 


-7 


-10 


-6 


-3 


-6 


-1 


-10 


-7 


-7 


-7 


-6 


-1 


-5 





-9 


-8 


-3 


-6 


-8 


-6 


-8 


-5 


-10 





1-4 


-3 


-5 


-6 


-6 


-10 





-6 


-9/ 
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The critical graph of this matrix consists of two s.c.c. comprising 6 and 3 nodes 
tively. They are shown in Figures H] and [5l together with their cychc classes. 

Figure 4. Critical graph of ([M]) 






/// 



IV 




Figure 5. Cyclic classes of the critical graph 



The components of C{A) induce block decomposition 



(63) 
where 



A 



All Ai2 
A21 A22 



(64) 



All 



/-8 -1 -8 
-4 -5 -2 
-7 -9 -8 



-8 -9\ 

-6 

-8 -4 

10 -7 -4 



7-4-8 
2 -7 -10 -6/ 



\0 -1 

The core matrix and its Kleene star are equal to 



A 



22 



(65) 



A 



Core 



Core\* 



-1 
-1 



/-5 

-5 



-9\ 

-10 
-6 -9 J 
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By calculating A, A^, ... we obtain that the powers of A become periodic after T{A) = 6. 
In the block decomposition of A^ analogous to (l63l) . we have the following circulants: 



4(6) 



(66) 



A. 



(6) 
21 



/ 


-1 


-2 





-1 


-2\ 


-2 





-1 


-2 





-1 


— 1 


—2 





— 1 


—2 








-1 


-2 





-1 


-2 


-2 





-1 


-2 





-1 


V-1 


-2 





-1 


-2 


0/ 




-1 


-2 


-3 


-1 


-2\ 


-2 


-3 


-1 


-2 


-3 


-1 


V-1 


-2 


-3 


-1 


-2 


-3/ 



(6) 
12 



A, 



(6) 
22 



The corresponding blocks of "reduced" power 

i{6) 

are 



A 



(6) 
11 



-2\ 



(67) 



A 



(6) 
21 






-1 


4(6) _ 
) ^12 — 


-1 


-2 


-1 


-2 








-1 


-2/ 


-1 


-2^ 






-3 


-2\ 


-3 


-1 


4(6) _ 

; ^22 — 







-3 


-2 


-V 




1-3 


-2 


0/ 



/-2 


-1 


-1\ 


-1 


-2 


-1 


-1 


-1 


-2 


-2 


-1 


-1 


-1 


-2 


-1 




-1 


-2/ 


U 


Q 

— 


-2\ 


-2 





-3 




-2 




-1 - 







Note that and are Kleene stars, with all off-diagonal entries negative. 
Using fl66|) . we specialize system (jH]) to our case, we see that this system of equations 
for the attraction cone Attr(A, 1) consists of two chains of equations, namely 



(6J 



Xi®X4® (Xg - 1) © (Xg - 1) = 

= X2 © X5 © (X7 - 1) © (Xg - 1) = X3 © X6 © (xj - 1) © (xg - 1), 
(X2 - 1) © (X5 - 1) ©X7 = 

= (X3 - 1) © (X6 - 1) © X8 = (Xi - 1) © (X4 - 1) © Xg. 



Note that only and —1, the coefficients of (^A'"°^^)* (which is equal to A'""'^'^ in our 
example), appear in this system. 
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7.3. Algorithm for the strongly connected case. Here we consider a 6 x 6 max-plus 
example 



/-3 





-1 


-19 


-7 


-3\ 


-3 


-4 





-10 


-19 


-16 





-3 


-1 


-10 


-8 


-8 


-1 


-4 


-4 


-1 


-1 


-3 


-1 


-1 


-4 


-2 


-4 


-1 


V-4 


-2 


-4 


-1 


-4 


-1/ 



and apply to it the algorithm described in Subsect. 16.31 The critical graph of this matrix 
consists just of one cycle of length 3, and there are 3 non-critical nodes. 



2 




1 < 3 > 



Figure 6. Critical graph and non-critical nodes of (!69l) 
The core matrix in this case is equal to 



/o 


-10 


-7 


-3\ 


-1 


-1 


-1 


-3 


-1 


-2 


-4 


-1 


\-2 


-1 


-4 


-V 



Vector h = (—10 — 7 — 3)"^, whose components are computed by 

3 

(70) ^i = ^ ^ki, for i = 4, 5, 6, 

k=l 

comprises 2, 3, 4- components of the first row of A'~^°^'^. The arguments of maxima in fITOl) 
give, after the cyclic shift by one position, the Boolean vectors 

(71) P4 = (1 1), P5 = (0 1 0), Pe = (0 1 0). 

These vectors encode, for the corresponding non-critical nodes t = 4,5, 6, the starting 
cyclic classes (here, just critical nodes!) of paths which go from C{A) directly to t and 
whose length is 3. 
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B 



/-I -1 -3\ 

-2 -4 -1 
[-1 -4 -l) 



B* 



/ -1 -2\ 

-2 -1 

v-1 -2 o; 



The lengths of optimal non-critical paths (whose weights are entries of B*) can be written 
in the matrix of sets 



(72) 



U 



/{O} {1} {2}\ 
{1,2} {0} {1} 

V{1} {2} {0}/ 



Further we compute 



0B* = (-10 -7-3) 



/ -1 -2\ 

-2 -1 
[-1 -2 o) 



(-4 -5-3) 



The maxima in 0^ htb^^ for all i are achieved only at t = 6, so W4 = W5 = Wq = {6}. 
Hence G4, G5 and Gq are shifted Pq and the shift is determined by the components in the 
last row of U which is ({1} {2} {0}). From Pq = {0 1 0) we conclude that 



G4 = (0 1), G5 = (1 0), Ge = (0 1 0). 



By Proposition 16.81 we have that Gi{t) = 1 if and only if t G K^^\i) (where r > T{A) 
is a multiple of 7 = 3). Using this rule we obtain that K^'^\l) = {5}, K^^\2) = {6}, 
f^('")(3) = {4}^ and using the vector of coefficients B* = (—4 — 5 — 3), we can write 
out the system for attraction cone 



(73) 



Xi © (Xs - 5) = X2 © {Xq - 3) = X3 © (X4 - 4). 
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To verify this result, we observe that in our case T{A) = 8 and 
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-3 
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-1 
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-2 


-2 


-1 


-5 


-7 


-5 


-1 
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-6 
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1-2 


-2 


-3 


-7 


-7 





Applying cancellation to the critical subsystem of ^ x = ® x, we obtain (1731) . 
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